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Based on mean field calculations with Skyrme interactions, we extract a constraint on the isovector 
effective mass in nuclear matter at saturation density po , i.e., m*(po) = (0.77d=0.03)?n by combining 
the experimental data of the centroid energy of the isovector giant dipole resonance (IVGDR) and the 
electric dipole polarizability «d in 208 Pb. Meanwhile, the isoscalar effective mass at po is determined 
to be m*(po) = (0.91 d= 0.05)ra by analyzing the measured excitation energy of the isoscalar giant 
quadrupole resonance (ISGQR) in 208 Pb. From the constrained m*(po) and m*(po), we obtain the 
isospin splitting of nucleon effective mass in asymmetric nuclear matter of isospin asymmetry S 
at po as [m*(po,d) — m*(po, <5)]/m = Amt(po)S + 0((5 3 ) with the linear isospin splitting coefficient 
Ami(po) = 0.33=b0.16. We notice that using the recently corrected data on the cm in 208 Pb with the 
contribution of the quasideuteron effect subtracted slightly enhances the isovector effective mass to 
ra*(po) = (0.80=b0.03)m and reduces the linear isospin splitting coefficient to Ara*(po) = 0.27±0.15. 
Furthermore, the constraints on m*(p), m*(p) and Ara*(p) at other densities are obtained from the 
similar analyses and we find that the A m\ (p) increases with the density. 

PACS numbers: 21.65.-f, 24.30.Cz, 21.60.Jz, 21.30.Fe 


I. INTRODUCTION 

Nucleon effective mass, which is usually introduced to 
characterize the dynamical properties for the propagation 
of (quasi)nucleons in nuclear medium, is of fundamental 
importance in nuclear many-body physics dm. While 
there exist several different kinds of nucleon effective 
masses in non-relativistic and relativistic approaches [H- 
0, we shall focus in this work on the total nucleon effec¬ 
tive mass used typically in the non-relativistic approach, 
which measures the momentum dependence (or equiva¬ 
lently energy dependence by assuming an on-shell disper¬ 
sion relation) of the nucleon single-particle potential in 
nuclear medium. In isospin asymmetric nuclear matter, 
neutrons and protons may feel different single-particle 
potentials which can then lead to the isospin splitting 
of nucleon effective mass, i.e., m*_ p = (m* — m*)/m. 
The isospin splitting of nucleon effective mass may have 
a profound impact on various physical phenomena and 
quantities in nuclear physics, astrophysics and cosmol¬ 
ogy [§, 9], such as the properties of mirror nuclei 101, 
transport properties of asymmetric nuclear matter [11 — 
0 , neutrino emission in neutron stars 0, and the pri¬ 
mordial nucleosynthesis in the early universe [22]. The 
isospin splitting of nucleon effective mass is also related 
to the momentum dependence of the nuclear isovector 
(symmetry) potential in nuclear medium [23] and thus 
the nuclear symmetry energy [2~l] [28j[ which is of critical 
importance for many issues of both nuclear physics and 
astrophysics but remains largely uncertain. A most re¬ 
cent review on the isospin splitting of nucleon effective 
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mass as well as its relation to the symmetry energy and 
symmetry potential can be found in Re. 0. 

Theoretical studies based on either microscopic many- 
body theories or phenomenological approaches have thus 
far given widely divergent predictions on m*_ p . For ex¬ 
ample, non-relativistic Brueckner-Hartree-Fock and rela¬ 
tivistic Dirac-Brueckner-Hartree-Fock calculations indi¬ 
cate m*_ p > 0 [29M3l| in neutron-rich matter, while 
relativistic mean field, Skyrme-Hartree-Fock (SHF) and 
Gogny-Hartree-Fock models predict either m* _ p > 0 or 
ra* < 0 [ 6 . Ehi l32l - [37l | , depending on the interactions. 
During the last several years, significant progress has 
been made in determining the isospin splitting of nu¬ 
cleon effective mass by analyzing experimental data [9]. 
However, there is still no quantitatively and even qualita¬ 
tively consensus on the behavior of m* _ p in asymmetric 
nuclear matter. For example, while the optical model 
analyses of nucleon-nucleus scattering data [24], [381 favor 
m*_ p > 0 in neutron-rich matter at po, the transport 
model analysis on the double n/p ratio in heavy ion col¬ 
lisions seems to suggest the opposite conclusion [39j (but 
see Ref. [ 13 ]). Therefore, any new and independent con¬ 
straints on the isospin splitting of nucleon effective mass 
are extremely helpful for understanding the issue on the 
behavior of m*_ p in asymmetric nuclear matter. 

Nuclear giant resonances provide an important ap¬ 
proach to determine nucleon effective mass. It has been 
well established that the excitation energy E x of the 
isoscalar giant quadrupole resonance (ISGQR) in finite 
nuclei is related to the isoscalar effective mass m* (p) (nu¬ 
cleon effective mass in symmetric nuclear matter) at po, 
i.e., m* )0 (see, e.g., Refs. jul-Eij). A value of m* 0 ~ 0.8m 
has been estimated by analyzing experimental data for 
ISGQR excitation energy in early studies [Hj, and more 
recent microscopic random phase approximation (RPA) 
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calculations suggest that the ISGQR in heavy nuclei fa¬ 
vors m* o ^ 0.9m HM. Moreover, within the RPA 
approach using Skyrme interactions, the isovector effec¬ 
tive mass m*(p) (i.e., neutron (proton) effective mass 
in pure proton (neutron) matter) at po, he., m* 0 , is 
closely related to the enhancement factor k of the en¬ 
ergy weighted sum rule (EWSR) mi in the isovector gi¬ 
ant dipole resonance (IVGDR) [451-147] . Unfortunately, 
while the peak of the IVGDR strength function has been 
well located for a number of nuclei by photoabsorbtion 
measurements j48|, neither the mi nor the k has been ac¬ 
curately determined. Furthermore, the detailed relation 
between m* 0 and mi or hi has not yet systematically in¬ 
vestigated for different nuclei. Therefore, the m* 0 has so 
far not yet been properly constrained. 

The properties of the heavy doubly magic nucleus 
208 Pb, especially including various kinds of giant reso¬ 
nances, have been well researched. In the present work, 
we mainly study how the ISGQR and IVGDR of 208 Pb 
constrain the m*_ p in neutron-rich matter. Thanks to 
the recent high resolution measurement for the electric 
dipole polarizability a d in 208 Pb, which is determined by 
the inverse energy weighted sum rule m_i of the IVGDR, 
at the Research Center for Nuclear Physics (RCNP) (49j . 
in this work, we deduce the mi of the IVGDR in 208 Pb 
from the experimental value of the IVGDR centroid en¬ 
ergy E -1 = i/mi/m-i [48|. Using the RPA calculations 
with a number of representative Skyrme interactions, we 
establish the detailed relations between m* 0 (m* 0 ) and 
the E x of ISGQR (mi) in 208 Pb, and then extract rel¬ 
atively accurate constraints on m* and m* from the IS¬ 
GQR excitation energy E x and the mi of the IVGDR 
in 208 Pb, respectively. Within the SHF model, we show 
that the m*_ p is completely determined by the m* and 
m*, and thus we can obtain constraints on the m*_ p , 
which is the main motivation of the present work. For 
the first time, our results indicate that the data on the 
giant resonances in 208 Pb definitely favor m*_ p > 0 in 
neutron-rich matter, which would be very helpful to pin 
down the isospin splitting of nucleon effective mass. 


II. MODEL AND METHOD 

A. Nucleon effective mass in Skyrme-Hartree-Fock 
approach 


the neutron/proton Fermi momentum, U q is the single¬ 
nucleon potential, and e q is the nucleon single-particle 
energy satisfying the following dispersion relation 

£q{k,p,S) = E- + U q (k,e q (k),p,5). (2) 

In this work, we use the standard Skyrme interaction 
with a zero-range and velocity-dependent form as [47] 


Ui 2 (R,r) = £ 0 (1 +x 0 P a )S(r) 

+ ];t 3 (l + x 3 P a )p°(R)6(r) 

6 

+ ^(l + x^CK'H^+S^K 2 ) 

+ t 2 (l + x 2 P a )K .8(r)K 
+ iW 0 (a! + a 2 ) ■ [K x<S(r)K], (3) 


with r = ri — r 2 and R = ( 1*1 + r 2 )/2. In the above 
expression, the relative momentum operators K = (Vi — 
V 2 ) /2i and K' = — (V 1 — V 2 )/2 i act on the wave function 
on the right and left, respectively. The quantities P a and 
Gi denote, respectively, the spin exchange operator and 
Pauli spin matrices. In the following, several Sky rme 
interactions with nonstandard spin-orbit term [50| are 
also employed, but the spin-orbit term is irrelevant to 
the expressions of nucleon effective mass. 

Within the standard SHF approach, the nucleon effec¬ 
tive mass in asymmetric nuclear matter with density p 
and isospin asymmetry 5 can be expressed as [47j 


h 2 1 

2m* (p, <S) 2 m 4 1 
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1 + 0*1 P - o + *1 )Pq 


1 


1 

+ -t 2 


1 + -X 2 P+ - + X 2 pq 
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(4) 


By setting p q = p/2 in Eq. (JU), the isoscalar effective 
mass can then be obtained as [47]] 


h 2 

2 m* 3 (p) 


2m 


+ + — ^(4x 2 + 5)p. 


(5) 


The isovector effective mass, which corresponds to the 
proton (neutron) effective mass in pure neutron (proton) 
matter, can be obtained with p q = 0 in Eq. (j4]) as [13 


In non-relativistic approaches, the effective mass m* of 
a nucleon q (n or p) in asymmetric nuclear matter with 
density p and isospin asymmetry 5 = (p n — p p )/(p p + p n ) 
can be calculated as [1] 


m q {p , S) 
m q 


1 + 


m q dU q (k , e q (k, p, 6),p, S) 


dk 


kl 


, (1) 


where m q represents the mass of neutrons or protons 
in free-space (m q = m is assumed in this work), k q F is 


2^) - i + ^ l(xi + 2)p+ l h{X2 + 2)p • (6) 

From Eqs. dU), (|5j) and (|6]), one can obtain the isospin 
splitting of nucleon effective mass, i.e., 
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where the isospin splitting coefficients A^ 2 n _i(p) can be 
expressed as 


Am* 2n _ 1 (p)=2 


ml 



(8) 


The above expressions reveal that, within the SHF 
model, the ra*_ p is completely determined by the m* 
and m*, and the sign of ra*_ p in neutron-rich matter is 
the same as that of m* — m*. 


B. Random-phase approximation and nuclear giant 
resonances 


The random-phase approximation 0 provides a suc¬ 
cessful microscopic approach to study giant resonance ob¬ 
servables in finite nuclei. Within the framework of RPA 
theory, for a given excitation operator Fjm 5 the reduced 
transition probability from RPA ground state |0) to RPA 
excitation state \v) is given by: 


B(EJ : 6 —> \v)) = |<i / ||-Fj||6)| 2 


2 


J2(X^ + Y^)\(m\\Fj\\i) 


(9) 


where m(i) denotes the unoccupied (occupied) single nu¬ 
cleon state; (m\\Fj\\i) is the reduced matrix element of 
Fjm'i and and Y^ are the RPA amplitudes. The 
strength function then can be calculated as: 


S(E) = J2MFj\\ 0)\ 2 S(E-E u ), (10) 

V 


where E v is the energy of RPA excitation state | v ). Thus 
the moments of strength function can be obtained as: 

m k = f dEE k S(E ) = ]T |(i,||Fj||6>| 2 £t (11) 

J 7/ 


For the IVGDR and IVGQR that we are interested in 
here, the excitation operators are defined as: 


N Z 

Am = ^Xk <y i M ( fi ) - J Y riYlM (^)’ ( 12 ) 

i= 1 i= 1 

A 

F 2M = J2 r i Y 2M(fi), (13) 

i=t 

where Z, N and A are proton, neutron and mass number, 
respectively; is the nucleon’s radial coordinate; Iim(^) 
and Y2M.(fi) are the corresponding spherical harmonic 
function. 


C. Nucleon effective mass and nuclear giant 
resonances 

It is well known that the isoscalar effective mass at 
saturation density, i.e., m* 0 , is intimately related to the 


excitation energy of the ISGQR in finite nuclei. In the 
harmonic oscillator model, the ISGQR energy is IS IH 


F x 



(14) 


where ftwo is the frequency of the harmonic oscillator. 
This semiempirical expression reveals the correlation be¬ 
tween the ISGQR excitation energy and the isoscalar ef¬ 
fective mass m* 0 , which has been also confirmed by mi¬ 
croscopic calculations [H, Q . 

Meanwhile, the isovector effective mass at saturation 
density m* 0 is correlated with the energy weighted sum 
rule mi of the IVGDR [52|, i.e. 


m i 


9 h 2 NZ 
47r 2 m A 


(1 + k). 


(15) 


where k is the enhancement factor reflecting the devia¬ 
tion from the Thomas-Reiche-Kuhn sum rule [53| (e.g., 
due to the exchange and momentum dependent force). 
Within the Skyrme-RPA approach, k is given by [13, 




WWz / 

' h ( 1 + t) + * 2 ( 1 + f). • 


(16) 


Substituting Eqs. (|6]) and (fl6j) into Eq. (fl5l) leads to 


mi 


9 h 2 NZ 
47r 2 m A 
A 


1 + 


NZ 



\ /pn(r)pp(r)rf 3 r 

) Po 


(17) 


which suggests that the EWSR mi (and thus k) of the 
IVGDR is proportional to for a fixed nu¬ 

cleus. In particular, by assuming p n = p p = po/2, one 
then obtains the following approximate expressions [47| 


9_tf_NZ_ /m^r 1 
4w 2m A \ m J ’ 


and 


1 /( 1 + k). 


(18) 


(19) 


III. RESULTS AND DISCUSSIONS 


To study the correlation between the nucleon effec¬ 
tive mass and the giant resonance observables, we se¬ 
lect 50 representative Skyrme interactions [H], EH EH 
(i.e., BSkl, BSk2, BSk5, BSk6, BSkl3, Es, Gs, KDE, 
KDEOvl, MSk7, MSLO, MSL1, NRAPR, Rs, SAMi, SGI, 
SGII, SK255, SK272, SKa, SkI3, SkM, SkMP, SkM*, 
SkP, SkSl, SkS2, SkS3, SkS4, SkSC15, SkT7, SkT8, 
SkT9, SKX, SKXce, SKXm, Skxsl5, Skxs20, SLy4, 
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SLy5, SLylO, SV-K241, v070, v075, v080, v090, vl05, 
vllO, Zs, Zs*). The corresponding ISGQR excitation en¬ 
ergies and EWSRs of the IVGDR in 208 Pb are calculated 
by using the Skyrme-RPA program by Colo et al [52]. 

In the calculation of the ISGQR excitation energy E x , 
we smear out the strength function with Lorentzian func¬ 
tions with a width 1 MeV. We note that varying the 
width has little influence on the peak energy. The ob¬ 
tained data-to-data relations between 10 3 /F 2 in 208 Pb 
and m s ^/m predicted by the chosen 50 Skyrme interac¬ 
tions are displayed in Fig.HJ Also included in Fig.[l]is the 
linear fit together with the corresponding Pearson corre¬ 
lation coefficient r. As expected from the semiempirical 
relation Eq. (HI, one can see that a strong linear correla¬ 
tion exists between 1/F 2 and m* 0 /m with the coefficient 
r as large as 0.971. And the linear fit gives 


10 3 

= (0.66 ± 0.26) ± (8.49 ± 0.30) 
Et 



( 20 ) 


where the E x is in MeV. 

In the present work, we invoke the weighted average 
of experimental values for the ISGQR energy in 208 Pb, 
i.e., E x = 10.9 ±0.1 MeV [43j] , which is shown as the 
hatched band in Fig. [TJ Combining this w< ' 
and Eq. m, we extract the isoscalar ef 
saturation density as 


—^ =0.91 ±0.05. 
m 

Here the error is obtained from the proj 

experimental uncertainty of E x and para 

the linear fit. This constraint is consiste: 

suit m* o ~ 0.8 — 0.9m obtained from an 

GQR of Nd and Sm isotopes [56[, and nat 

the empirical value of m* 0 ^ 0.9m prec 

Skyrme interactions which are obtained ffi 

perimental data of the ISGQR excitation 

nuclei [44],|45|. It is also in good agreement 

of m* o ~ 0.92m from the extended Brueckner-Hartree- 

Fock calculation with realistic nucleonic forces [29]. 

For the IVGDR, we use the chosen 50 Skyrme inter¬ 
actions to evaluate the mi of the IVGDR in 208 Pb with 
energy up to 130 MeV. Similarly, in Fig. [2j we plot the 
data-to-data relations between 10 4 /mi and m* 0 /m as 
well as the linear fit and Pearson correlation coefficient 
r. It is clearly shown that an excellent linear correlation 
exists between 1/mi and m* 0 /m, and the linear fit gives 


10 4 r yx\> ^ 

-= (2.17 ±0.05) ±(11.5 ±0.07)—(22) 

mi m 


where the mi is in MeV-fm 2 . Experimentally, the cen¬ 
troid energy of the IVGDR, i.e., E-\ = \Jrrii/m-i, in 
208 Pb has been well determined from photoabsorbtion 
measurements, i.e., E_ 1 = 13.46 MeV j48[, and the in¬ 
verse energy weighted sum rule m_i can be obtained 
from the experimental value of the electric dipole polariz¬ 
ability measured at RCNP, i.e., a d = 20.1 ± 0.6 fm 3 . 



FIG. 1: (Color online) 10 3 /F 2 in 208 Pb vs m*^/m predicted 
by a large number (50) of Skyrme interactions. The linear fit 
gives 10 3 /F 2 = (0.66±0.26) + (8.49±0.30) ( m*^/m ) with the 
Pearson correlation coefficient being 0.971. The hatched band 
corresponds the weighted averages of the experimental values 
for the ISGQR excitation energy in 208 Pb, E x = 10.9 ±0.1 
MeV 0. 
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FIG. 2: (Color online) 10 4 /mi in 208 Pb vs m%^/m predicted 
by a large number (50) of Skyrme interactions. The linear 
fit gives 10 4 /mi = (2.17 ± 0.05) + (11.5 ± 0 . 07 ) 771 *^ 0 /^ with 
the Pearson correlation coefficient being 0.999. The hatched 
band (cyan band) corresponds to the (corrected) experimental 
value of the EWSR mi of the IVGDR in 208 Pb (see text for 
the details). 


through the following simple relation 

m -' = (23) 

The experimental values of E- 1 and m_i together thus 
give mi = 905.60 ± 27.03 MeV • fm 2 . Therefore, one can 
constrain isovector effective mass at saturation density 
using Eq. ([22]) . and the result is 

m* o 

—^ = 0.77 ±0.03. 

rn 


( 24 ) 
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One can see that our constraint is rather accurate and 
well consistent with the empirical value, e.g., m* Q /m = 
0.90 ± 0.2 from analyses of finite nuclei mass data (57j |. 
In addition, from the well known relation Eq. (fl5j) for 
the EWSR, the value of the enhancement factor k can 
be deduced as k = 0.228 ± 0.037 with m± = 905.60 =b 
27.03 MeV • fm 2 , which is in very good agreement with 
k = 0.22 dz 0.04 reported in Ref. [58| and consistent with 
the estimate of k « 0.2 — 0.3 in Ref. We note that 
using the relation m* 0 /m ~ 1/(1 ± k) (i.e., Eq. (Il9]) i 
leads to a little bit larger ^asK~ 0.30 ±0.05, indicating 
that Eq. (HU) is indeed satisfied approximately. However, 
one should be cautious to use the relation m* 0 /m ~ 
1/(1 ± k) for an accurate determination on m* 0 from /c, 
and vice verse. 

From the constraints on m* 0 and m* 0 , one can 
then obtain the isospin splitting ra*_ p (po) according to 
Eq. (|5]). In particular, we obtain the first-order (linear) 
isospin splitting coefficient Am*(p) at po as 

Am*(p 0 ) =0.33 ±0.16, (25) 

which is in very good agreement with the constraint 
Am*(po) = 0-32 ± 0.15 obtained in Ref. (24| and the 
more recent constraint Ara*(po) = 0.41 ± 0.15 extracted 
in Ref. [38[ from the global optical model analysis of 
nucleon-nucleus scattering data. The present result is 
also consistent with the Ara*(po) = 0.27 obtained by an¬ 
alyzing various constraints on the magnitude and den¬ 
sity slope of the symmetry energy at po poj. The 
positive value of Araf(po) further agrees with the mi¬ 
croscopic Brueckner calculations with realistic nuclear 
forces 1291431). In addition, it is interesting to see that 
the higher-order isospin splitting coefficients are rather 
small and can be neglected safely. For example, the third- 
order isospin splitting coefficient Ara^po) is found to be 
0.01 ± 0 . 01 . 

The above analyses are only made at saturation den¬ 
sity po and it is also interesting to see the constraints on 
m*(p), m*(p) and m*_ p (p) at other densities. Similar 
analyses indicate that the strong linear correlation also 
exists between 1 /E% and m*(p)/m as well as between 
1/rai and m*(p)/m at other densities p. Shown in Fig. [3] 
(a) are the constraints on the m*(p)/m and m*(p)/m 
as functions of density extracted from the ISGQR and 
IVGDR in 208 Pb, respectively. The inset of Fig. [3] (a) 
shows the density dependence of the corresponding Pear¬ 
son correlation coefficient r for 1/F 2 vs ra* (p )/m as well 
as 1/rrii vs m*(p)/m. The corresponding constraints on 
the isospin splitting coefficients Am*(p) and A m^(p) as 
functions of density are shown in Fig. [3] (b). Also in¬ 
cluded in Fig. [3] (b) are the Ara*(po) constraints ob¬ 
tained in Refs. [24], [38|, [6Q| as discussed earlier. Indeed, 
one can see that all the r values in the inset of Fig. [3] 
(a) are larger than 0.95 for 0 < p < 0.32 fm -3 that we 
are considering here, indicating the strong linear corre¬ 
lation. Particularly, the strongest correlation appears at 
p « 0.19 fm -3 (with r = 0.97035) for 1/F 2 vs m*(p)/m 
while at p ^ 0.13 fm -3 (with r = 0.99961) for 1/rai vs 
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FIG. 3: (Color online) Panel (a): Constraints on the density 
dependence of the isoscalar and isovector effective mass, ra* 
and ra*, extracted from the ISGQR and IVGDR in 208 Pb, 
respectively. The inset shows the corresponding Pearson cor¬ 
relation coefficient r as a function of density. Panel (b): 
Constraints on the density dependence of the isospin split¬ 
ting coefficients Amf(p) and A rris(p) obtained in this work. 
The hatched bands (cyan bands) represent the results of ra*, 
Ara*(p) and A rris(p) without (with) subtracting the contribu¬ 
tion from the quasideuteron effect. The AraJ(po) constraints 
obtained in Refs. HI EE El are also included for comparison. 


m*(p)/m. It is seen that the m*(p)/m is generally larger 
than m*(p)/m and both m*(p)/m and m*(p)/m decease 
with density but the latter exhibits a stronger density de¬ 
pendence, which leads to the isospin splitting coefficients 
Am*(p) and A m^(p) increases with density as observed 
in Fig. [3] (b). It is interesting to see that the third-order 
isospin splitting coefficient A m\(p) is very small (about 
0.05 even at p = 0.32 fm -3 ) and can be approximately 
negligible. The stronger isospin splitting of nucleon ef¬ 
fective mass at higher densities may have implications 
on the isospin effects in heavy ion collisions and neutrino 
emission in neutron stars as mentioned earlier, and these 
are deserved further explorations in future. 

Very recently, Roca-Maza et al. pi] pointed out that 
the measured value of ckd in 208 Pb reported in Ref. |49] 
is contaminated by the non-resonant quasideuteron ef¬ 
fect at higher energies. The quasideuteron effect should 
be mainly due to correlated neutron-proton pairs in nu¬ 
cleus and it overwhelms the IVGDR for the photon ab¬ 
sorption at higher energies beyond about 25 MeV [62l |. 
In principle, the contribution from the quasideuteron ef- 
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feet should be subtracted to directly compare the ex¬ 
perimental strength against the theoretical RPA cal¬ 
culations for the IVGDR. In Ref. [61] (and references 
therein), this contribution has been determined and the 
electric dipole polarizability in 208 Pb has been corrected 
to be 19.6 zb 0.6 fm 3 . Invoking the corrected value of 
c^d, we derive a value of the EWSR of the IVGDR as 
mi = 883.05 =b 27.03 MeV • fm 2 which is plotted as the 
cyan band in Fig. [2j Repeating the above analyses, one 
can further obtain the isovector effective mass at po as 
m* 0 /m = 0.80 ± 0.03 and the corrected linear isospin 
splitting coefficient Am*(p) as Am*(p 0 ) = 0.27 =b 0.15. 
The obtained new value of the enhancement factor is 
k = 0.197 db 0.037, which is still in good agreement with 
the results in Refs. jH, [59|. Again, we extract the den¬ 
sity dependence of the isovector effective mass m* and 
the isospin splitting coefficients Am*(p) and Am 3 (p) by 
using the corrected value of an in 208 Pb due to the 
quasideuteron effect, and the results are shown as the 
cyan bands in Fig. [3j Overall, one can see that the cor¬ 
rection on the experimental value of an in 208 Pb from 
subtracting the contribution of the quasideuteron effect 
does not affect the qualitative conclusions and only leads 
to small corrections on quantitative results. 

Furthermore, we have made similar analyses for the 
IVGDR of the semidouble-closed-shell nucleus 68 Ni. Us¬ 
ing the measured centroid energy E-\ — 17.1 ± 0.2 
MeV [63] as well as the electric dipole polarizability 
an = 3.88 ± 0.31 fm 3 obtained in Ref. [ 6 l| from a 
Lorenzian(-plus-Gaussian) extrapolation of the measured 
GDR strength [63] to the high-energy (low-energy) re¬ 
gion, we extract an isovector effective mass at po as 
m * o /m = 0.81 ± 0 . 11 , which is in good agreement with 
the above results obtained from analyzing the data of 
208 Pb although the uncertainty is larger. In addition, 
we would like to point out that the present RPA cal¬ 
culations are not expected to reproduce the experimen¬ 
tal spreading width of the GDR, and this problem can 
be solved effectively by taking into account the coupling 
to the collective low-lying (mainly surface) vibrations or 
phonons [64l-[67|. As discussed in Ref. [ 6 l|, however, such 
an effect beyond the mean-field approximation is not ex¬ 
pected to significantly affect the integral properties of the 
calculated strength that we are focusing on here. 

IV. CONCLUSIONS 

Based on mean field calculations with Skyrme inter¬ 
actions, we have demonstrated that the isoscalar and 
isoveetor effective masses at saturation density, i.e., m* 0 
and ra* 0 > can be well constrained by the ISGQR ex¬ 
citation energy E x and the EWSR mi of the IVGDR 


in 208 Pb, respectively. In particular, invoking the ex¬ 
perimental data for E x in 208 Pb, we have obtained the 
constraint m* 0 = 0.91 ± 0.05m. Meanwhile, combining 
the experimental IVGDR centroid energy and the electric 
dipole polarizability an = 20.1 ±0.6 in 208 Pb, we have de¬ 
duced a value of mi = 905.60 ± 27.03 MeV • fm 2 , and fur¬ 
ther extracted a value of m* 0 = 0.77 ± 0.03m. From the 
extracted m* 0 and m* 0 , we have obtained a constraint 
on the first-order (linear) isospin splitting coefficient of 
nucleon effective mass, i.e., Am*(po) = 0.33 ±0.16 which 
is in good agreement with the constraints extracted from 
global nucleon optical potentials constrained by world 
data on nucleon-nucleus scattering [24, |38| and is also 
consistent with the value obtained by analyzing the con¬ 
straints on the symmetry energy |60| . 

Furthermore, we have constrained the isoscalar and 
isoveetor effective masses as well as the isospin splitting 
of nucleon effective mass at other densities by the similar 
analyses of the giant resonances in 208 Pb. Our results in¬ 
dicate that the isospin splitting of nucleon effective mass 
increases with the density, and the third-order or higher- 
order isospin splitting coefficients are negligibly small. 

In addition, we have also investigated how our results 
change if the recently corrected experimental value of 
an = 19.6 ± 0.6 fm 3 in 208 Pb due to the quasideuteron 
effect is used. Our results indicate that the corrected 
value leads to mi = 883.1 ± 27.0 MeV • fm 2 , m * 0 = 
0.80 ± 0.03 m, and Am*(po) = 0.27 ± 0.15. Therefore, 
the quasideuteron effect in 208 Pb only plays a minor role 
on the extractions of the isoveetor effective mass and the 
isospin splitting coefficient of nucleon effective mass. Our 
present work reveals for the first time that the data on 
the giant resonances in 208 Pb definitely favor m* > m* 
in neutron-rich matter, which sheds a light upon under¬ 
standing the isospin splitting of nucleon effective mass in 
asymmetric nuclear matter. 
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